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Abstract. We present a new method for computation of the 
Korteweg-de Vries hierarchy via heat invariants of the 1-dimensional 
Schrodinger operator. As a result new exphcit formulas for the 
KdV hierarchy are obtained. Our method is based on an as- 
ymptotic expansion of resolvent kernels of elliptic operators due 
to S. Agmon and Y. Kannai. 



1. Introduction and main results 

1.1. Heat asymptotics of Schrodinger operator and the KdV 
hierarchy. Consider the 1-dimensional Schrodinger (or Sturm-Liouville) 
operator 

Its heat kernel H{t, x, y) is a fundamental solution of the heat equation 

The heat kernel has the following asymptotics on the diagonal as t — > 0+ 
(see [3]): 

I oo 

H{t,x,x)r^^=y2hr,[u]e, (1.1.2) 

where hn[u] are some polynomials in u{x) and its derivatives. The coef- 
ficients hn[u] are called heat invariants of the 1-dimensional Schrodinger 
operator. 

Computation of heat invariants of self-adjoint elliptic operators is a 
well-known problem in spectral theory which has many applications, 
in particular to geometry and theoretical physics (see [3, 4, 6]). Heat 
asymptotics of the 1-dimensional Schrodinger operator are of particular 
interest due to their relation to the Korteweg-de Vries (KdV) hierarchy 
which is one of the basic objects in the theory of integrable systems 
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(see [9, 10]). Namely, the KdV hierarchy is defined by ([2]): 



I - I'^-M. (1.1.3) 



where 



(2nV 

Gn[u] = jr^hnH neN. (1.1.4) 

Set Mo = u, Un = d'^u/dx"', n G N, where Un, n > are formal vari- 
ables. The sequence of polynomials G„[m] = Gn[uo,Ui,U2, ■ ■ ■] starts 
with ([2]): 

Gi[u] = uo, 6*2 [m] = U2 + 3ul, Gslu] = M4 + lOuoU2 + 5ul + IOmq, . . . 

(1.1.5) 

In particular, substituting 6*2 [^] into ( |1.1.3|) we obtain the familiar 
Korteweg-de Vries equation ([9, 10, 5, 2]): 

du d^u g du 
dt dx^ dx 

1.2. Explicit formulas for the KdV hierarchy and the heat in- 
variants. There are many ways to compute the polynomials Gni^] 
and /in[w]. Originally it was done through a recursive relation between 
the heat invariants (see [5, 8]). However, in order to understand the 
structure of heat invariants and the KdV hierarchy closed expressions 
are desirable. This might also make actual computations easier. Such 
expressions were first obtained by R. Schimming and I. Avramidi ([13, 
2]), all of them having a very complicated combinatorial structure. 

In this paper we present a method for computation of heat invariants 
of the 1-dimensional Schrodinger operator (see [11] for a similar method 
for the Laplacian on a 2-dimensional Riemannian manifold). As a 
result we obtain simpler explicit formulas for /i„[m] and hence for the 
KdV hierarchy. Our method is based on an asymptotic expansion of 
resolvent kernels of elliptic operators due to S. Agmon and Y. Kannai 
([1, 11])- 

We introduce the following notation. Consider L^[x'^^ ^"]|^=q, k > 
n>l. This is a polynomial in m(0),m'(0), . . . , ^'•^^"^-'(O). Let us make a 
formal change of variables: u^^\0) Ui, i = 0, . . . ,2k — 2, and denote 
the obtained polynomial by Pkn[u]- 

Theorem 1.2.1. The polynomials Gn[u], n G N, are given by: 



2n m 



Gn[u] = E E 1 Ltl):. p^M. (ra 



^n.=„22-2n+i(2™)(2A;-2n)! 
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Expanding the expression for (see [12]) we get: 
Theorem 1.2.2. The polynomials Gn[u\, n G N are equal to: 

2n m 



m=n k=n \ 2n Ji,---,jp 

^ ' Jl + '--+jp = 2(n-p) 

0<«o</i<---<ip-i=fc-p VJl/ \ -'2 / \ JP 

2/i>jiH hji+i, i=0,...,p-l. 

Expressions for the heat invariants /i„[u] can be easily deduced from 
([TXTI ) and (|i;23) using (|1X§. 



Theorems |1.2.1| and |1.2.2| are proved in sections 3.1 and 3.2. 



2. Computation of heat invariants 

2.1. A modification of Agmon-Kannai expansion. The original 
Agmon-Kannai theorem ([1]) deals with asymptotic behaviour of re- 
solvent kernels of elliptic operators. In [11] we have obtained a concise 
reformulation of this theorem suitable for computation of heat invari- 
ants. Before presenting it we introduce some notations. 

Let if be a a self-adjoint elliptic differential operator of order s 
on a Riemannian manifold M of dimension d < s and let Hq be the 
operator obtained by freezing the coefficients of the principal part H' 
of the operator H at some point x G M: Hq = H'{x). Denote by 
Rx{x,y) the kernel of the resolvent R\ = {H — A)~^, and by Fx{x,y) 
— the kernel of Fa = {Hq - \y\ 

Theorem 2.1.1. ([11]). The resolvent kernel Rx{x,y) has the follow- 
ing asymptotic representation on the diagonal as X ^ oo: 

oo 

Rx{x,x) ~ p(x)-i E XmF;^+\x,x), (2.1.2) 

m=0 

where p{x)dx is the volume form on the manifold M and the operators 
Xm are defined by: 

= pyi)'(^^H'H^-', m > 0. (2.1.3) 
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2.2. Prom the resolvent to the heat kernel. Let us return to our 
particular case of 1-dimensional Schrodinger operator. We have 



M = W 



d 



H 



-L, Hn 



p{x) 



1. 



Note that ( p.l.2|) is an expansion in powers of —A due to the following 
formula ([1]): 



dx" 



Fr\x,x) = {-xy 



2lT 



m+1 



(2.2.1) 



Collecting terms with the same powers of —A in ( 2.1.2 ) one receives 
the standard asymptotic expansion of the kernel of Rx = {—L — X)~^ 
on the diagonal as A — oo (see [1]): 



R,{x,x)^i~X)-'/'Y.^n[u]{-X) 



-nil 



(2.2.2) 



n=0 



where r„[M] are polynomials in u[x) and its derivatives. 

For each ra G N the polynomials r„[M] and the heat invariants /?.„[«] 
differ by a constant multiplier (cf. [11]): 

(2n)! 



-Mn\v\ 



22n+l^! 

3. Proofs of main theorems 
3.1. Proof of Theorem [Qj] . Due to ( 11X41) and (g^lJ), 



(2.2.3) 



(3.1.1) 



Therefore in order to compute Gn[u] we need to find the coefficient 
r„[n] in (|X|). 

Without loss of generality set x = 0. We need to collect all terms 
in the sum ( 2.1.2|) ) containing (— A)^"/^. Let us apply formula ( 2.2.1 



Note that if a is odd then ( |2.2.1| ) is zero identically. For the terms we 
are interested in we have a = 2m — 2n and hence m > n. On the other 
hand the order of the differential operator Xm defined by ( p.l.3|) is not 
greater than m (see Lemma 5.1, [1]). Hence a = 2m — 2n < m which 
implies m < 2n. 

Denote by a polynomial in m(0), m'(0), m"(0), ... which is obtained 
from the polynomial rn[u\ by a formal change of variables Ui — u^^\0), 
i > 0. 

Then due to (|2.1.3| ) we have: 



2n m 
m=n k=n 



\rri+k 



^2m-2k 
^^2m,—2k 



X 



2m— 2n 



{2m-2n)\ 



x=0 
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where 

^mn = / 7-^^ T-TT = ' B im - n + -■,n + - 

271 i-oo (^2 + 271 V 2' 2 

and the integral above is computed using the tables ([7]). 
Let us rewrite the B-function in terms of factorials: 

1 1\ r (m-n + i) r (n + I) (2n)\ ■ (2m - 2n)\ ■ 7t 



2' 2 J r(m + l) 22"^ -m! -n! ■ (m-n)!' 

Substituting this into the expression for $mn and simplifying we 
finally get: 

where Pkn[u] = Pkn[uo, Ui, .., U2k-2] are polynomials obtained from L^[x'^^~'^"] 
by a formal change of variables 

m(0) Mo, m'(0) Ml, . . . , M^'^^-'HO) ^ M2fc-2. 

By ( |3.1.1| ), formula ( p.l.2| ) implies ( |1.2.1| ) and this completes the proof 
of the theorem. □ 

Remarks. The proof of Theorem |1.2.1| is similar to the proof of The- 
orem 1.4 in [11]. Formula ( |1.2.1|) was programmed using Mathematica 
([14]) and for 1 < n < 5 the results agreed with the already known 
ones (cf. ( [TX5| ), [5]). 

3.2. Proof of Theorem |1.2.2| . In [12] an explicit formula for the 
powers of the 1-dimensional Schrodinger operator was obtained. We 
may rewrite it in the form: 

^2k k ^jo 

= x:^ + E E C',i,...,,,M,,(x)..M,^(x)— -, 

p=l Jo>h'---'ip "-^ 

iO+JlH hjp=2{*:-p) 

where 

C. . = J2 /2/oU2/i-jiy _ /2/p-i-ji jp-i\ 

0<lo<h<-<lp-i=k~-p \Jl/\J2j\ Jp J 

2li>ji-\ hji+i, i=0,...,p-l. 

Setting jo = 2k— 2n this allows to expand the expression for [x2^~2n.j 
and hence for the polynomials Pkn[u]- Substituting it into (|3.1.2| ) we im- 
mediately receive ( |1.2.2| ) which completes the proof of the theorem. □ 
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